Electronic states and magnetic excitations in LiV204: 
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Abstract. Motivated by recent inelastic neutron scattering experiment we 
examine magnetic properties of LiV2 04. We consider a model which describes 
the half-filled localized A\ g spins interacting via frustrated antiferromagnetic 
Heisenberg exchange and coupled by local Hund's interaction with the 1/8-filled 
itinerant E g band, and study it within an exact diagonalization scheme. In the 
present study we limited the analysis to the case of the cluster of two isolated 
tetrahedrons. We obtained that both the ground state structure and low-lying 
excitations depend strongly on the value of the Hund's coupling which favors 
the triplet states. With increasing temperature the triplet states become more 
and more populated which results in the formation of non-zero residual magnetic 
moment. We present the temperature dependence of calculated magnetic moment 
and of the spin-spin correlation functions at different values of Hund's coupling 
and compare them with the experimental results. 



PACS numbers: 75.25,+z, 75.30.Et, 75.30.Vn, 71.45.Lr 
1. Introduction 

In recent years, there has been a great interest in the study of the origin of the heavy 
fermion (HF) behavior observed in the paramagnetic transition- metal oxide L1V2O4. 
This compound has a spinel structure and its pyrochlore lattice consists of corner- 
shared tetrahedrons of V +3 5 ions located in slightly distorted oxygen octahedron. It 
is the first metal showing heavy fermion behavior without any / orbital pQ. The 
electronic specific heat coefficient 7 = C e /T = 0.42 J/Mol K 2 at T = 1.0 K is the 
highest value measured for 3d transition-metal oxides. 

In addition to having quasi-particles with an unusually heavy effective mass, 
L1V2O4 has some peculiar magnetic properties. The magnetic susceptibility and 
inelastic neutron scattering measurements indicate a spin liquid behavior over a 
large range of intermediate temperatures E]. The absence of magnetic order 
indicates that geometrical frustration, due to the pyrochore structure, is relevant in 
this system 4 . On the other hand, the recent neutron scattering experiments [5] reveal 
that in addition to antiferromagnetic correlations, ferromagnetic-like correlations on 
V sites appear over some temperature range. These observations are inconsistent 
with the 'classical' spin-liquid picture, and demonstrate that itinerant contribution of 
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strongly correlated electrons, which could lead to the effective ferromagnetic exchange 
interaction between localized spins is also important. 

There were some suggestions to explain HF behavior in LiV204 using analogy 
with 4/— systems where these effects are attributed to the hybridization of localized 
4/— levels and itinerant spd bands. This analogy has been based on the fact that 
according to the band structure calculations Ej the trigonal distortion splits t2 g 
orbitals into singlet A\ g and doublet E g states, and the center of A\ g band lies lower 
than that of E g band. Therefore, one can make a plausible assumption and treat 
A\ g level as being occupied by a localized electron, whereas consider E g doublet as a 
quarter filled conducting band 0. 

In this paper we study the magnetic properties of L1V2O4 within the exact 
diagonalization analysis of a small cluster. The model that we consider includes 
both purely Heisenberg-like contribution from the super-exchange interaction among 
localized spins and an effective ferromagnetic double-exchange contribution driven by 
the itinerant electronic excitations. 

In order to capture some charge fluctuation of the system, we solve the model 
exactly for the cluster consisting of two disconnected tetrahedrons with either 5 and 
7 electrons, or only with 6 electrons, projecting out higher energy states with 4 and 8 
electrons. 

2. The Model 

In order to consider some charge fluctuations, we study a cluster of two isolated 
tetrahedrons, i.e. we don't take into account either the super-exchange interaction 
between localized spins of different tetrahedrons or the hopping of electrons between 
them. For this cluster, we can write the total Hamiltonian as H = Hi + H2, where 
-ffi(2) is defined on a single tetrahedron: 

-J„^S«< + j£S t a S«. (1) 

Here a — I, 2 is the tetrahedron's index; are ^-spins representing the localized Ai g 
electrons; c°J Q { c j a a) ^ ne creation (annihilation) operator of an itinerant electron with 
spin a =|, I and orbital a = 1, 2, corresponding to the ^-doublet. 
The first term in the Hamiltonian Eq. (JIJ describes the electron hopping between 
the nearest neighbor V ions, t"? being the transfer amplitude. The second term 
concerns the Hund's coupling Jh between the localized spins and the local spin 
density er,; of the itinerant ^-electrons. Finally, the last term describes the nearest- 
neighbor antiferromagnetic super-exchange interaction J among localized spins. Here 
we consider infinite on-site Coulomb repulsion between itinerant electrons and project 
out states with double occupancy. 

Considering only direct overlap between 3d wave functions \(dda) — —0.281 eV, 
(ddn) — 0.0076 eV] |7], transfer matrix elements tf? can be easily calculated using 
the Table I of the paper by Slater and Koster 

In further analysis, we consider the value of Hciscnbcrg exchange coupling J — 10 meV 
which is in agreement with the estimates given in the literature ^2 El The Hund's 
coupling Jh is known to be little screened in solids, and can be simply related to its 
atomic value. The estimates reported in the literature suggests that Jh = 0.68 — 1.0 
eV [H El El El We use in our study the value J H = 0.8 eV obtained by the LDA+U 
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ab initio calculations [Sj, however we also consider the variation of Jh over a wider 
range. 

2.1. Hilbert space sectors 

To define the Hilbert space, we choose a basis in which each state |n) is a product 

\n) = \N 2> S 2> a 2 h®\N 1 ,S 1 ,a 1 ) 1 , (2) 

where the state ■ • -) a characterizes the tetrahedron a — 1,2. The index N a is the 
number of electrons (localized plus itinerant) for the tetrahedron a, with an averaged 
charge of 6 electrons per tetrahedron N% + N 2 = 12. For a given tetrahedron, the spin 
state <S a = X)i=i characterizes the total spin of the four localized electrons. These 
four spins can be coupled either to singlet, triplet or quintet states. If Hund coupling 
is neglected, considering only the antiferromagnetic super-exchange interaction, the 
ground state consists of two-fold degenerate singlet with energy — 3J, followed by 
ninefold degenerate triplet states and five quintet states with the energy 3J. The 
quantum number a a describes all the remaining degrees of freedom, in our case, 
the spin state and the orbital and site occupancy of the itinerant electrons on the 
considered tetrahedron. 

In order to calculate both the ground state and the finite-temperature properties 
of the cluster of two isolated tetrahedrons, first, we performed a numerical exact 
diagonalization of the Hamiltonian Hi characterizing one tetrahedron. Because of 
large degeneracy, the Hilbert space of one tetrahedron consists of 256, 1536 and 
4096 states for 5, 6 and 7 electrons, respectively. However, since the Hamiltonian 
is invariant under rotation of the total spin, it is also diagonal by blocks in the basis 
of the eigenstates of the total spin (itinerant plus localized). 

2.2. Probabilities 

For a given temperature T = 1//3, the average value of an operator S (e.g. the nearest 
neighbors spin-spin correlation) is defined by the usual relation 

(5} = i2<n|ea:K-/M0S|n). (3) 

n 

Choosing a basis diagonalizing cither the operator S or the Hamiltonian, this relation 
can be rewritten as 

(S) = T)<n|E|n>. (4) 

n 

Here, the density matrix p(n,T) = -g- (n\exp(~ (3H)\n) is the probability for the state 
\n) to be occupied at temperature T. Using the definition Eq. @ of the state \n) of 
the cluster of two isolated tetrahedrons, the probability P(N, S) that the tetrahedron 
1 occupies a state with N electrons and a local spin state S is: 

P(N,S)=^ Yl {N,S, ai \{\2- N,S 2 ,a 2 \exp{-PH)\\2- N,S 2 , a2 )\N,S, ai ) 

Ql ,02,52 

= ^z?LY J {N,S, ai \exp{-pH 1 )\N,S,a 1 ), (5) 
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Figure 1. a) Ground state spin probabilities at different values of Hund's 
coupling Jff. b) Charge probabilities at temperatures T = 0.1,0.3,1.0 eV for 
J n = 0.8 eV. c) and d) Spin probabilities at temperatures T = 2, 10, 50 meV for 
Jn = 0.0 eV and Jjj = 0.8 eV, respectively. 



where the partition function for one isolated tetrahedron with TV = 5, 6 or 7 electrons 
is defined as Zjq = J2s a (N,S,a\exp(—(3Hi)\N,S,a). The partition function of the 
two tetrahedrons system can be cast into 

Z = Z 5 Z 7 + Z 6 Z 6 + Z 7 Z 5 . (6) 

For the further analysis, it is convenient also to define the probability that the 
tetrahedron 1 occupies a state with a total local spin S: 

p s MS) = J2p(n,s), (7) 

N 

and the probability to occupy a state with TV electrons: 

P charg e(N) = J2 p (N,S). (8) 
s 



3. Results and discussion. 



3.1. Low energy states 

First, let us discuss the weight of the different spin subsectors in the ground state. 
Fig-Da) illustrates the spin probabilities P(N = 6,5) for T = and different values 
of Hund's coupling. At Jh — 0.0 eV (decoupled spins and conducting electrons), 
the ground state is a two-fold degenerate singlet. As soon as the Hund's coupling 
is switched on, the spin degeneracy of the ground state is lifted and some triplet 
components occur. The weight of the triplet states is more than 20% for a realistic 
coupling Jjj = 0.8 eV, but the singlet components remain important even for much 
larger coupling Jjj — 2.0 eV. For the three values of coupling that we considered, 
the quintet contributions to the ground state are negligible. In Fig. Q] a) we present 
only the probabilities of states corresponding to the tetrahedron with TV = 6 electrons 
because the states with 5 and 7 electrons are higher in energy and, therefore, all the 
probabilities P(5,S) and P(7,«S) are equal to zero. At low temperatures, as in the 
ground state, mainly the states with TV — 6 electrons are occupied. The charge sectors 
with TV = 5 and N — 7 contribute only at rather high temperatures. In Fig.^b) the 



Electronic states and magnetic excitations in LiV204.' Exact diagonalization study. 5 




T eV T eV 



Figure 2. The temperature dependence of a) the specific heat C e and b) the 
coefficient 7 = C e /T. On both a) and b), the temperature axis is logarithmic. 

temperature evolutions of P c harge{N) are presented for Jh = 0.8 eV. However, we 
should note that only static charge fluctuations are taken into account here, and 
considering interacting tetrahedrons would certainly increase the charge fluctuations 
contribution at low temperatures. 

The evolution of the spin probabilities P S pin{S) with increasing temperature are 
presented in Fig. ^ c) and d) for Jh = 0.0 eV and Jh — 0.8 eV, respectively. 
The weights of the triplet states become substantial at T ~ 10 meV, and the spin 
probabilities are mainly proportional to the spin degeneracies at T ~ 50 meV. 

To conclude the description of the low energy excitations, let us discuss the 
thermodynamical properties such as the electronic specific heat C e and the coefficient 
7 = C e /T. In Fig. [3 a) and b) we present the temperature dependence of C e and 
7, respectively. For all values of Hund's coupling, the specific heat has three peaks 
at finite temperatures. Two low-energy peaks correspond to spin excitations. For 
a non zero Hund's coupling, the energies of the magnetic states are decreased, and 
the corresponding peak occurs at a lower temperature. The third and the highest 
temperature peak characterizes the charge sector excitations and it is independent of 
the Hund's coupling. 

3.2. Local moment and magnetic correlation junctions 

We discuss now the temperature dependence of the spin correlation functions at 
different values of Hund's coupling. We denote as S'°* = Si+ai the total spin at a given 
V ion. In Fig. a), we plot the nearest-neighbor spin correlation function (5'°' S^J) ■ 
At low temperatures the correlations are antiferromagnetic for the three values of Jh 
considered. When the temperature is increased the correlations change from antiferro- 
to ferro-type at a temperature T ~ 30 meV for Jjj = 0.8 eV and Jh = 2.0 eV, 
whereas it is always antiferromagnetic for Jh — 0.0 eV. To emphasize the role of the 
itinerant electrons, we plot on the inset of Fig. 03 a) the nearest-neighbor correlation 
functions (SizSW) of the localized spins only, which remain antiferromagnetic at any 
temperature and Jh- 

This observation is in qualitative agreement with the inelastic neutron scattering 
measurements of Murani et al. 5: where the correlations are antiferromagnetic (with 
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Figure 3. a) Temperature dependence of the total spin nearest neighbors 
correlation functions. The inset shows, for comparison, the correlations of only 
the localized spins, b) Temperature dependence of the local magnetic moment. 
On both a) and b), the temperature axis is logarithmic. 



a wave vector around Q = 0.6 A -1 ) below T ~ 2 K, and develop ferromagnetic- 
like correlations (corresponding to a peak at Q = 0) with increasing temperature. 
A possible interpretation of these data could be the following: at low temperatures 
the ferromagnetic exchange, induced by the itinerant electrons due to the double 
exchange mechanism |16| . is weaker than the direct antiferromagnetic exchange, and 
the resulting correlations are antiferromagnetic. With increasing temperature, carriers 
become more mobile, and, as a consequence, the effective ferromagnetic exchange 
grows. In Fig. [3] b) we present the local effective moment per site, which is defined 
as: m e ff(m e ff + 1) = ((XL S*°*) 2 )/4. At Jh — 0.0 eV, when the localized spins are 
completely decoupled from the itinerant electrons, the magnetic moment is determined 
only by the localized spins. At zero temperature, the ground state is a singlet, and 
magnetic moment is zero. With increasing temperature low lying triplet states become 
thermally populated, and this causes the formation of the local moment. At high 
temperature the obtained value of 0.67 corresponds to a paramagnetic moment of 1.5 
electrons in average on a V ion. 

At Jh = 0.8 eV and Jh = 2.0 eV, the ground state has some triplet components but 
the weight of the singlet is still important. As a consequence, the effective moment is 
small at low temperature. 

Finally, we calculated the static magnetic susceptibility, defined as xCO = 
((Sf 0t ) 2 )/T. At high-temperature, for all values of the Hund's coupling the inverse 
susceptibilities showing linear behavior , corresponding to a Curie- Weiss law. The 
Curie- Weiss temperatures extrapolated from high temperatures (T ~ 500 — 1000 K ) 
are 9 = -551 K, -475 K and -360 K for J H = eV, 0.8 eV and 2.0 eV, respectively. 
A detailed analysis of the magnetic susceptibility measurements Q] gives 9cw ~ —37 
K when the fit is performed over the range T ~ 100 — 300 K, but it gives 6 ~ —600 
K if fitted over the range T ~ 500 — 1000 K. This latter value is of the same order as 
the one we have calculated, and we guess that the corresponding range of temperature 
is such that the correlations between tetrahedrons are negligible. According to our 
calculations, this temperature does mainly characterize the spin correlations but in 
the charge sectors N — 5 and N = 7. The value dew ~ —37 K, fitted at lower 
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temperatures, would be more characteristic of spin correlations in the N = 6 charge 
sector. 

4. Conclusion 

We have performed an exact diagonalization of the small cluster, consisting of two 
disconnected tetrahedrons. The size of the cluster, constraint on the number of elec- 
trons and no exchange between two tetrahedrons gives us the possibility to make some 
qualitative description of the system, its spectrum and magnetic behavior. 
Wc find that the ground state is a two fold degenerate singlet when the contribution 
from the itinerant electrons is not considered. The degeneracy is lifted and some im- 
portant triplet components appear when the itinerant electrons are coupled to the 
localized spins by the Hund's exchange interaction with a realistic value. Consistent 
with experimental observations [5] , the obtained temperature dependence of the spin 
correlation function shows a crossover from antiferromagnetic to ferromagnetic be- 
havior when the temperature is increased. The increasing population of the magnetic 
states with temperature results in the formation of a non-zero residual magnetic mo- 
ment which is also observed in the experiments. 

We thank A. P. Murani, J. R. Iglesias and G. Jackeli for valuable discussions. 
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